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Magnetic Polarizability of a Short Right Circular 

Conducting Cylinder l 

T. T. Taylor 2 

(August 1, 1960) 

The magnetic polarizability tensor of a short right circular conducting cylinder is calcu- 
lated in the principal axes system with a uniform quasi-static but nonpenetrating applied 
field. One of the two distinct tensor components is derived from results already obtained in 
connection with the electric polarizability of short conducting cylinders. The other is calcu- 
lated to an accuracy of four to five significant figures for cylinders with radius to half-length 
ratios of %, }&, 1, 2, and 4. These results, when combined with the corresponding results for 
the electric polarizability, are applicable to the problem of calculating scattering from cylin- 
ders and to the design of artificial dispersive media. 

1. Introduction 

This article considers the problem of finding the magnetic polarizability tensor fi i5 of a 
short, that is, noninfinite in length, right circular conducting cylinder under the assumption 
of negligible field penetration. The latter condition can be realized easily with available con- 
ducting materials and fields which, although time varying, have wavelengths long compared 
with cylinder dimensions and can therefore be regarded as quasi-static. The approach of the 
present article is parallel to that employed earlier [1] 3 by the author for finding the electric 
polarizability of short conducting cylinders. 

The magnetic polarizability tensor relates the induced magnetic dipole moment m t to the 
inducing field Bj according to the equation 

mi=PijB j} (1) 

with the cylinder situated in free space. The centered coordinate system of principal axes 
used in [1], in which the :-axis coincides with the axis of rotational symmetry, is also employed 
here. In this system, the equal transverse components fi xx and fi vy are denoted /3 tt , and the longi- 
tudinal component fi zz is denoted fin. In agreement with other conventions of [1], a represents 
the radius of the cylinder; 6, the half-length; and c, the half -diagonal or distance from center to 
edge, nationalized M.K.S. units are used throughout. 

The method described has been used successfully for calculating fi tt for cylinders with a/b 
ratios of }{, %, 1, 2, and 4. For the calculation of fin, the following relation, proved in appendix 
A, is invoked: 

/•oAi— 5^ (2) 

where a tt is the transverse electric polarizability treated in [1]. It is interesting that eq. (2) 
is valid not only for circular cylinders but for all conducting solids of revolution. Both fin 
and fi tt are negative, indicating that an artificial medium composed of short cylinder elements 
arranged in a lattice will be diamagnetic. 

In the transverse case, where no simple relationship such as (2) is known to exist, the prob- 
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lem which presents itself is that of finding, in the external space, a static magnetic field which 
has no normal component at the surface of the cylinder and which reduces to a uniform field at 
infinity. With such a field there is associated physically a surface current density j on the cylin- 
der and a field free condition in the interior. This surface current is related to the field values 
at the surface by : 

M oj = e w XB, (3) 

where e n is the outward normal, and it gives rise to the induced magnetic dipole moment : 

m=±jrX)ds' (4) 

As with the corresponding problem involving electric polariz ability, there is no known coordi- 
nate system in which a short cylinder can be treated by the conventional method of separating 
Laplace's equation and a different method, involving arbitrarily good approximations to the 
current distribution function, must be employed. 

In a manner analogous to that of [1], let the cylinder be regarded merely as a geometrical 
construct to which the surface current distribution is firmly affixed. If the uniform (applied) 
component of the magnetic field is subtracted, the remaining field will be that due to the cur- 
rent system on the cylinder acting alone ; it will be characterized by a dipole-like external field 
and a perfectly uniform internal field which is equal and opposite to the applied field. The prob- 
lem may therefore be restated in terms of finding, on a finite cylindrical surface with closed 
plane ends, a surface current distribution which, when acting alone, generates a uniform mag- 
netic field in the interior. The technique employed here makes use of this viewpoint and con- 
sists in setting up, on the surface, a current density function which is completely determined 
by a finite number of parameters, then solving for the values of these parameters so as to obtain 
maximal uniformity of the current field or, alternatively, maximal cancellation of the applied 
field, within. 

The research reported here, like that related to the electric polarizability, was performed 
under the direction of Prof. W. K. Smythe. The author is indebted to him for many valuable 
discussions and to the personnel of the California Institute of Technology Computing Center 
for their cooperation. 

2. The Transverse Problem 

In this problem, let a uniform magnetic field of magnitude B be applied in the positive 
^/-direction. The induced surface current density j may be regarded as equal to J£ where £ is 
the thickness of an infinitesmal surface layer in which the current is confined and J is the 
value of the conventional current density vector averaged over the thickness of this layer. 
In the quasi-static case, J is solenoidal and can be regarded as the curl of an equivalent aver- 
aged magnetization M which is also confined to the infinitesmal surface layer and which is 
normal to the surface at all points. On the side of the cylinder, let 

if^=-aF(s/&)sin*, (5) 

and on the upper end, 

i^l = - a ^(p/ a ) s in^>, (6) 

where F and G are dimensionless functions as yet to be determined. The configuration just 
described constitutes a nonuniform magnetic shell, that is a magnetic shell [2, p. 261] associated 
with a distributed current system. 
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Evidently the surface currents on the side are given by 



^=F(z/b) cos *, 



(7) 



H)J* 



f- s =~F'(z/b)sm<l>, 



(8) 



where the prime denotes differentiation with respect to the dimensionless argument z/b. The 
surface currents on the upper end are 



Ms = _«<jf(p/ O )cos0; 
n p 



(9) 



^=fl"G,/a)sintf. 



(10) 



The current distribution on the lower end is equal and opposite to that on the upper end. 
Typical forms for the functions involved in eqs (7) through (10) are illustrated in figure 1; a 
sketch of the actual lines of current flow is presented in figure 2. 
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(Q/p)G(p/Q) 



-* p/a 



G'(/o/a) 




■► p/Q 




Figuke 1 . Typical forms for the functions describing the equiv- 
alent magnetic shell and the current system. 



Figuke 2. Sketch of the current system. 



The asymptotic behavior of current and field near the edge of the cylinder is determined 
by the fact that, locally, the edge resembles a configuration consisting of two infinite perfectly 
conducting planes intersecting at right angles with the field restricted to the outer (larger) of 
the two angular domains. Such a configuration is a familiar one in two-dimensional potential 
theory. In the case at hand, two mutually perpendicular components of magnetic field are 
distinguished: (a) the component which is everywhere parallel to the edge, and (b) the com- 
ponent which bends around the edge. Each of these field components has associated with it 
a surface current obeying eq (3). The component of magnetic field parallel to the edge is 
locally uniform in the outer space and is of no further interest in the present discussion. The 
surface value of the component which bends around the edge, however, increases as the edge 
is approached and is in fact proportional to Z" 1/3 where I is the perpendicular distance from 
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the edge to the surface point in question. This dependence, which is the same as that of the 
electrostatic charge density on a similar freely charged configuration, is obtained by means of 
the Schwarz-Christoffel transformation. The associated component of surface current is, of 
course, parallel to the edge and everywhere proportional to the field. On the cylinder, then, 
the (^-components of current, and hence the functions F r and G f must be asymptotic to con- 
stants times l~ lfZ as I tends to zero. This, in turn, imposes certain requirements upon F and 
6. Quantitatively, 

^(^/6)-ft+^ 2 (6-0) 2/3 as z^b; (11) 



G(p/a)~Cz-C 4 (a-p) 2 



as 



>a; 



(12) 



~F'{zlb)~^C 2 (b-zY 



1/3 



as 



(13) 



G'(p/a)~-giC A (a—py 



1/3 



as 



(14) 



It is seen that two edge conditions, C z = Ci and C±=C 2) must be satisfied. The former of these 
insures continuity of current flow across the edge of the cylinder, the latter that the <j> com- 
ponents on side and end match one another as both tend to infinity at the edge. 

The functions F and G are represented by finite weighted sums involving the unknown 
coefficients f b , j m , g b , and g m , and the \{/ m functions defined in [1]. Note that the parity para- 
meters are chosen so as to make F even and G odd and that the parameter which controls the 
singularity is given the value 2/3. 
Thus, 



F(z/b) = * 



and 



G( P /a) = 



The constant terms 



o, 



a m=o \ 6 a/ 



A=/^o(^o,o,|), 



M<&| 

p<a> 
P>a. 



(15) 



(16) 



(17) 



P „7 



9vz=9»f<{ J'i'O 



(18) 



are called the "basic" terms as in [1], Together they constitute a very simple current system 
in which the flow takes place in plane rectangular loops around the cylinder, the plane of each 
loop being parallel to the x^-plane of figure 2. Superposed upon this basic flow pattern is 
another pattern described by the nonbasic (summed) terms in (15) and (16). The addition 
of this pattern gives a gentle curvature to the current lines and introduces the necessary singular 
behavior at the edges. Clearly the first edge condition demands that g b =f b and henceforth 
the single symbol f b will be used as the coefficient of the basic term in both F and G. 

In relating the self field of the current system to M£, a magnetic scalar potential U, defined 
such that B=— VU, will be used. The boundary conditions at the magnetic shell are the 
duals of the corresponding boundary conditions for the electric potential at a charged surface, 
that is, the normal gradients on the two sides of the shell are equal but the potential values 
differ by the local intensity ol the magnetization. Quantitatively, if I and II are two adjoining 
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regions separated by a magnetic shell, and if e n is a unit vector pointing from I into II, then 

[/il-Dl=/zoM w f. (20) 

Application of these boundary conditions and of other well-known procedures of potential 
theory leads to the following expression for the potential in the interior of the cylinder due to 
the current system on the side: 

U 8 (p,4> 9 z) —a 



cB 



=-^ r R f F(z'/b) cos kz'dz'l kaK[{ka)I x {kp) cos kz dk sin </>. (21) 



The Fourier cosine transform of the f function, enclosed in the square brackets above, is easily 
expressed in terms of Bessel functions with the aid of eq (A9) in [1]. The function h{kp) cos 
kz sin <fi, which is regular in a neighborhood of the origin, may be represented in a spherical 
harmonic expansion about the origin. Thus (21) becomes 

U,(rM) ^ ~ab (2\» C- Yr Jm(B) , *&=> , (-\) m J, m+m+ Akb) 



» -ab /2\" 2 ("" T J m (B) ^ (-l) m J tm+ m+Akb) -\ 

Pi c W J„ L'» (kby*^£='o Jm (kby*+' J 



cB p=o 

kaK[{ka) (kcy +i dk ^^ (!) 2?+lp k +>( cos V sin <t>. (22) 

Here and elsewhere, the singularity parameter v' for the nonbasic terms is understood to be 
equal to 2/3. The spherical harmonic expansion of (22) converges for all r<Ca. 

With reference to the ends of the cylinder, the interior potential generated by the current 
distributions on the two ends is given by 



U e (p,<l),z)_a 



cB 



fi ( a G(p'/a) J l {kp f )p'd P '~\ kae- kb J,{kp) cosh kzdk sin </>. (23) 



This, in turn, may also be expressed as an expansion in spherical harmonics which converges 
for all r<6: 



U e (r : 



,(r,fl,0) _^ £rV. Mka) , ^ n (-l) m J2m+2 + Aka) l 
cB ~^ e Jo b n (ka) "^ 0m (ha) 1 *' J 



1 / r \2p+l 

kae- kb (kcy^hlk {2p+2)l (y PL + i(cos 6) sin d>. (24) 



Evidently the sum of (22) and (24) constitutes an expansion of the interior potential due 
to the total current distribution, convergent within the largest sphere which can be inscribed 
in the cylinder. Such an expansion may also be expressed in the form 



Ns-l _ Ne-1 

cB p^o L m=o " 771^0 



TJ(r0d>) °° r iv ' 8 " 1 ^ e_1 "1 /r\ 2p+1 

-S I J* Z » + S f m Z7+f b Z?+ S ffw Z;-J (y PL + i(cos 0) sill </>, (25) 



where the Z's are obtained from (22) and (24) by integrating [3, p. 137] over k. These Z's 
become the matrix elements in a system of simultaneous equations which may be solved for 
the/ 6 ,/ m , and g m such that the coefficient of (rjc) 2p + l P\ v+l (cos 0) sin (/> will be unity for p=0 and 
zero for as many p^O as possible. The second edge condition, also included in the system ot 
simultaneous equations, is simply 

iVi— 1 //7\- 2 /3 Ne-l 

S /*+(?) S <7«=0. (26) 

m = W w=0 
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If one makes use of the fact that —kaK[(ka)=kaK (ka) J r Ki(ka), the expression for the 
matrix element Z s ™ becomes 

z ™=f (lT ( W+w£ [kaKo(ka)+Kl(ka)] J Tkiy^ ] W +ldh W 

Integration yields the following result : 

v.- (-l) m+ >V' +1 -"T(m+ P +i)T(m+p+i) a*/b\'»+' [Y , „ , 1\ 
r(|)r(2p+3)r(2m+|+/) bc\c) LV m ' tP ' t 2j 

2 F 1 (m+2>+|m- 2J +/;2m+|+/;|)+i 2 F l (m+i)+| m-p+l+v'; 2m+|+/;^)]- 

(28) 
The matrix element Z e p is simpler: 

z '--c w+w£ e ~ kt ^m^ 1 (/ " )2?+im - (29) 



This becomes 



7em _ (-l) m 2 2 " +1 -"T(m+p+i)V(m+p+2) b /a\ 2m+i 
r(i)r(2p+3)r(2m+34V) c\cj 

2 F 1 (m+p+^, m -p+l+v';2m+3+p';^- (30) 

The basic matrix elements, Z* 6 and Z e p b , are obtained by setting both m and v' equal to 
zero in the preceding formulas. 
Thus, 

^_ (-l) p 2 2?+1 r(ff+l)r(p+f) a 2 6[-/ l\ 

r(i)r(2i>+3)r(i) " <? \\ v+ 2) 

* (*+|-« I ^)4a^(p+|-l>+l; I ^)} (3D 
2^T(y+|)r(p+2) rf& / 5 a*\ 

z »" r(i)r(22>+3)r(3) c * 2 H^ + 2' p+ 1 > s ><*) (32) 

Using the appropriate formulas for analytic continuation [4, p. 108], one obtains 
et / , 5 , , a 2 \ r(3)r(-|) „/ , 5 , , 3 b 2 \ 

* {p+z-p+h 3; ^)= r( _;^; r (; + 2) °*v + 2'-* +i; 2 ; ?) 

1 e r(3)r(|) / 1 1.6«\ . 

+6r( ;P +f)r(-p+i)^ 1 (,-^+2'^+ 2 ' 2' ?/ (33) 



By [4, p. 3], 

T(-h) -(-l)T(p+|) 

r(-i»+i) r(f) 



(34) 



The second term in (33) vanishes unless p=0, in which case the hypergeometric function is 
easily summed : 

'■(KSK'-sr-e)- 

Let the symbol 8 0p be equal to unity for p=0 and to zero for p^O; then 

7«-» (-l)^+T(jt>+f)r( g +|) a^ / 5 _ , ..3.6^ . . 

z ° - So » r(i)r(2p+3)r(t) — "? »*V+2'-* +1 ' 2' ?/ (36) 
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In (25), the total coefficient of f b is the sum of Z 8 p and Z e p. This sum is denoted by Z b v and is 
given by 

Z "-^+ r(i)r(2p+3)r(t) " 7" L^ + 2> 2 H 2H ~2 ~^2' ?J 

4> (p+|-*+i ; I £)-£ (*+§)* (p+I'-p+I; | J)]- (37) 

Using the third Gauss relationship, as listed in [4, p. 103]: 

7»_, (-l)"2 2 T(p+f)r(y+l) a 2 6 „ / 3 3. 6 2 \ , v 

Z ""^ r(i)r(22»+3)r(f) 7 2jPi M ,+ 2 ,-2>+1 '2' ?/ (38) 

This has an even simpler form, obtained by an elementary transformation of the hypergeo- 
metric function and a simplification of the gamma function multiplier: 

The three formulas, (28), (30), and (39), contain the information necessary for a numerical 
calculation of the matrix elements. If both the electric and the magnetic polariz abilities are 
being calculated for the same cylinder, much labor can be saved by comparing the expressions 
for the matrix elements in this paper with those in [1] and exploiting the similarities. 
The magnetic dipole moment, given by eq (4), has only a 7/-component, 

m v=o ( — z 3*» sm *~~ a i« cos 4>)ad<t>dz+ b(j pe cos 4>— j 0e sin fycfopdp. (40) 

^J - b J JO Jo 

Using expressions (7) through (10), 

^ = ^f P* (z%F'sin 2 4>-aF cos 2 $)ad<t>dz + P P* b(~-G cos 2 4>-G' sin 2 4>) cfapdp; (41) 

=7ra 2 f I" f- (uF'-F)du- P (fl f +^fl f/ )^T (42) 



Moj^V 

5 



By parts, 



Ato my _ 

5 



i= T a 2 b\ [uF]-2 C Fdu-[uG] V (43) 

L o Jo J 

The first edge condition requires that F(1) = G(1); hence, 

^=-2Ta 2 b P Fdu. (44) 

The quantity m v jB is the transverse polarizability and the geometrical volume v is easily 
recognized. The orthogonality relationships given in [1] aid in the integration of F, yielding 
the following formula for the polarizability: 



Mo/3, 
^o 



"~^|/&+2 6/3 r( 18 ) *M ^^ 
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3. Kesults 

The two distinct tensor components of the magnetic polarizability are given in table 1. 
Actual values of the expansion coefficients for the F and G functions, which describe the 
surface current distribution via the equivalent magnetic shell for the transverse case, are 
included in appendix B. These data, upon which the f$ tt results of table 1 are based, were 
obtained by solving sets of 18 simultaneous equations. Experimentation with fewer equations 
has led the author to the conclusion that the accuracy is not quite as good as that obtained 
for the electric polarizability in [1]. In other words, the fl tt figures have a probable error of 
about 2 units in the least significant digit for a/b equal to unity and about 10 units in this 
digit for a/b equal to 4 or to }{. The fi u figures, on the other hand, were obtained from the 
electric polarizabilities simply by the application of eq (2), and therefore enjoy the same 
accuracy belonging to the latter. 

Table 1 



a/b 


Wfill/VO 


wfitt/vo 







-1.0000 


-2. 0000 


H 




-1.1575 


-1.8506 


v?, 




-1.3057 


-1.7351 


i 




-1.5853 


-1. 5795 


2 




-2. 1087 


-1.4140 


4 




-3.0907 


-1.2716 


CO 




— oo 


-1.0000 



v =2ira 2 b 



As an additional check upon the numerical calculations in the transverse problem, a 
procedure analogous to that employed in [1] was used. Consider the total field B, which is 
the sum of the applied field and the field due to the current system on the cylinder. A normal 
component B tn of the total field at any point on the surface of the cylinder indicates that an 
error is present. If n represents distance away from the surface in the outward normal direction 
at a particular surface point, one may calculate the displacement of the surface An which, con- 
sidering the local value of the normal field gradient, would reduce the normal field component 
to zero. In other words, let the following equation be solved for An: 



Btn" 



dn 



An=0. 



(46) 



This Aft may then be compared with the gross dimensions of the cylinder. The local field 
gradient bB tn fon depends upon the current distribution at the point in question and may be 
calculated by appealing to eq (3) and to the solenoidal character of B*. To first order, one 
has for any point on the side of the cylinder 



*B tp ( 



1 MnjL^hiX 

a d(t> dz J' 



(47) 



^ B sin <t> ( F J#_ F ,,\ 
a \ b 2 / 



&B tp B sin <f> 



(48) 



Let the p-component of the field due to the current system alone be denoted B„. Then 
B tp =B p -\-B sin <f>, An=Ap, and one obtains for the latter, 



B 



Ap B sin 4> 



fl 



-f+fi F" 



(49) 
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The quantities B p , F, F" and finally Ap/a were calculated for 2=0, that is for points on the 
"equator" of the cylinder; hence the term "equator check" in appendix B. 
The same method applied to the tipper end of the cylinder yields 



bB tz _ ( jt e bj^ 1 dj pe \ 

~te~-* ] \ P ~dp + p W 

5g fg= J3sin0 (_G^_ Q G\ 

bz a \ u u 2 / 

where u represents p/a. Substitution of the expansion AiU-\- A z u z ; + A 5 u 5 '+ 
is known to be an odd function, gives 



i>B t2 __ B sin <t> 
bz a 



(-%A z u— 24A 5 u z + . . .)• 



(50) 

(51) 
for G, which 

(52) 



Since both sides of (52) vanish at p=0, that is at the "pole" of the cylinder, one must differ- 
entiate with respect to p in order to obtain something useful. Thus, 



-*™L» (-8A 3 )= - iB t u V "(0). 



bpbz 



3a 2 



(53) 



Furthermore, for any point on the upper end, B tz is simply B z , the field due to the current 
system alone. Finally, 

Az 3a 2 b_ / B z 



b 4&G'"(0) dp VBsin0> 



(54 



Values of bB z /bp, G"', and Az/b were calculated for the pole; data are given under the heading 
of "pole check" in appendix B. Generally speaking, the pole checks are satisfactory but are 
not as good as the corresponding checks in [1]. 

4. Appendix A: Relation Between the Longitudinal Magnetic and the Trans- 
verse Electric Polarizabilities of a Conducting Solid of Revolution 

Let figure 3 illustrate a typical solid of revolution formed by rotating the curve p—p(z) 
about the 2-axis. 



s=s w is 

P \ ^^ 

/ \ 
L J ^ p 

s = 



Figure 3. Cross section of a typical solid of revolution 
showing coordinates and unit vectors. 



562759—60- 



207 



Let the unit vectors e s and e n point along the tangent and the outward normal of this curve, 
respectively. The variable s } associated with e s , measures distance along the curve and encom- 
passes the solid as it varies over the interval 0<s<S. It is sufficient if p(z) has a uniquely 
defined tangent and is normal at almost every point in this interval. 

If the applied fields consist of a uniform magnetic field B in the ^-direction and a uniform 
electric field E in the ^-direction, the following is to be proved : 



n rn z 



B 



1 Vx 

~2e E 



(Al) 



In the electric case, the scalar potential satisfies the conditions 

V=-Ef 1 ( P ,z) cos «p; 
/i=0 on the surface of the solid; 
/i~p at infinity; 
V 2 F=0. 



(A2) 



The magnetic case, on the other hand, is analyzed with the aid of the vector potential which, 
by symmetry, has only a ^-component. This component must vanish at s=0, s=S, and at all 
intermediate surface points; otherwise there would be a normal component of magnetic field 
at the surface. One finds that 



7? 

jf 2 =0 on the surface of the solid; 
/ 2 ~P at infinity; 
V 2 A=0. 



(A3) 



Application of the Laplacian operator to either cos </> ji or to e^/ 2 yields the same equation in p 
and z. In other words, 



{ 



_y 1 _&_!, &_ 

d P 2ir pd P P 2+ ds 2 



}{£}=»• 



(A4) 



Since j\ and f 2 obey the same differential equation, have the same boundary conditions at the 
surface of the solid, and have identical asymptotic behavior at infinity, it is clear that/i=/ 2 =/. 
The surface charge density in the electric case is 



<r=— e e n -vV=€ E cos <j> (e w .v/), 
and it is easily shown that the x-component of the induced electric dipole moment is 



p x =e Eir I (e n -vf)p 2 ds. 



The surface current in the magnetic case as given by (3) is: 

J=-e„X(vXA)=^e„X(vXeJ). 

Ho Zpo 

On the surface, where /=0, this reduces to the following: 
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(A5) 



(A6) 



(A7) 



J=-9-e*(ea-V/'), 



(A8) 



(A9) 



and the magnetic dipole moment is obviously 

m z = \ I (e n .yf)p 2 ds. 

Mo ^ Jo 
A comparison of (A9) with (A6) shows that the asserted relation (Al) is true. 

5. Appendix B: Values of the Expansion Coefficients for the Functions F and G 

Values of the expansion coefficients for the functions F and G and the results of the checking 
procedure at pole and equator are given in tables 2 through 6. The notation "A(p)" means 
"A times 10V 



Table 2 



a/b=H 



m 


fm 


Qm 


basic 


+0. 13640434 (+1) 

+0. 94335081 (+0) 
+0. 47094227 (+0) 
+0. 15962976 (+0) 
+0.52222212 (-1) 
+0.17272056 (-1) 

+0.55233089 (-2) 
+0. 16203994 (-2) 
+0.41694492 (-3) 
+0.90138919 (-4) 
+0. 15575188 (-4) 

+0.20048427 (-5) 
+0. 17022684 (-6) 
+0. 71319103 (-8) 






1 
2 
3 


-0. 37183474 (+0) 
-0. 18563192 (+0) 
-0.97767116 (-1) 


4 




5 




6 




7 




8 




9 




10 




11 




12 








Equator 
check 


Bp/E sin 0= 

Ap/a 


= -0. 99999999 
negligible 


Pole 

check. __ 


a-- (B z /B sin 0) = 

Op 

Az/b-- 


=+0. 11001902 
=+0.01246495 



Table 3 



a/b=H 



m 


fm 


Qm 


basic 


+0. 13390263 (+1) 

+0. 76790238 (+0) 
+0. 23079152 (+0) 
+0.58147254 (-1) 
+0. 16764875 (-1) 
+0. 50092196 (-2) 

+0. 13805064 (-2) 
+0. 32120205 (-3) 
+0. 57675668 (-4) 
+0. 70016001 (-5) 
+0. 42707647 (-6) 





1 

2 
3 

4 

5 
6 


-0. 32481666 (+0) 
-0. 14321008 (+0) 
-0.97128394 (-1) 
-0. 67455522 (-1) 
-0.36913235 (-1) 

-0. 11074155 (-1) 


7 




8 




9 








Equator 
check 


Bp/B sin 0—1. 00000004 
Ap/a : negligible 


Pole 
check. _ . 


a - - (B z /B sin 

Op 

A 


0) = -0.00614087 
2/6= -0.02256388 



Table 4 



o./6=l 



m 


fm 


Qm 


basic 


+0. 12992927 (+1) 

+0. 54323936 (+0) 
+0.92098198 (-1) 
+0.17619634 (-1) 
+0.38217146 (-2) 
+0. 67793221 (-3) 

+0.66769575 (-4) 





1 
2 
3 

4 

5 
6 


-0. 29402467 (+0) 
-0. 12410492 (+0) 
-0.80893285 (-1) 
-0.59953649 (-1) 
-0.44063915 (-1) 

-0.29238560 (-1) 
-0.16182197 (-1) 


7 




-0. 68464716 (-2) 


8 




-0. 19415309 (-2) 


9 




-0.27441144 (-3) 








Equator 
check. _ . 


B„/B sin 0= -0. 99999998 
Ap/a : neglible 


Pole 
check 


a ^- (Bz/B sin 0) = -0. 00000485 
Aijb= -0.00001512 



Table 5 



ajb =2 



m 


fm 


Qm 




+0. 12407197 (+1) 

+0. 33590275 (+0) 
+0.30150397 (-1) 
+0. 28485520 (-2) 






1 
2 
3 


-0. 25652350 (+0) 
-0. 11210208 (+0) 
-0.72413481 (-1) 
-0. 52386398 (-1) 


4 




-0.37934754 (-1) 


5 




-0.25696044 (-1) 


6 




-0. 15530579 (-1) 


7 




-0. 80562157 (-2) 


8 




-0. 34498633 (-2) 


9 




-0. 11632626 (-2) 


10 




-0. 28842547 (-3) 


11 




-0. 46632080 (-4) 


12 




-0. 36819245 (-5) 








Equator 
check. _. 




Bp/B sin 
ApJ 


£=-1.00000288 
a =+0.000001 15 


Pole 
check 


d 

a 5T 

dp 


(BJB sin 

Az/ 


) = -0.00000013 
6= -0.00000100 
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Table 6 



a/b=i 



m 


fm 


(Jm 


basic 


+0. 11775253 (+1) 

+0. 18232541 (+0) 
+0. 76983950 (-3) 






1 
2 


-0. 21277158 (+0) 
-0.98721317 (-1) 
-0.61346284 (-1) 


3 




-0.39649681 (-1) 


4 




-0.24392554 (-1) 


5 




-0.13622852 (-1) 


6 




-0.66941775 (-2) 


7 




-0.28211792 (-2) 


8 




-0.99369454 (-3) 


9 




-0.28363133 (-3) 


10 




-0. 62898095 (-4) 


11 




-0.10154428 (-4) 


12 




-0.10608227 (-5) 


13 




-0.53798542 (-7) 








Equator 
check 


Bp/B sin 0= -1.00378633 
A p /a = +0.00094822 


Pole 
check 


a — (B,/B sin </>) =+0.00000000 

Op 

Az/b: negligible 
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